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TWO-DIMENSIONAL ANALYTICAL AND EXPERIMENTAL PERFORMANCE COMPARISON 


FOR A COMPRESSOR STATOR SECTION WITH D-FACTOR OF 0.47 

by Nelson L. Sanger 
Lewis Research Center 

SUMMARY 

Analytically computed flow parameters for the midspan, double-circular-arc blade 
section of a highly loaded compressor stator have been compared with measured flow 
parameters. Analytical calculations were performed for inviscid flow on a blade-to.- 
blade flow plane, boundary layer on the blade suction and pressure surfaces, and total 
pressure loss coefficient. All velocities were subsonic. Comparisons were made at 
three incidence angles. 

Inviscid flow calculations require the outlet flow angle as input. Thus, either some 
trailing edge condition must be specified and the resulting exit angle accepted, or an exit 
angle may be specified and the resulting flow distribution accepted. For one method, 
iterations were performed on the exit angle until suction and pressure surface pressures 
were equal at the trailing edge. With this method, surface pressures diverged from 
measured values near the trailing edge, and exit angles differed by as much as 10° from 
measured values. A second method used measured values of exit angle as input to the 
inviscid flow calculations. Surface pressures showed much better agreement; pressures 
on suction and pressure surface became equal at between 91. 5 and 92. 5 percent chord 
for all incidence angles. 

Boundary layers were calculated as being turbulent from near the leading edge. 
Losses were calculated using a flow model developed by Speidel which includes a term to 
account for losses due to separation. To bring calculated losses into agreement with 
measured losses, very large boundary layer initial thicknesses were necessary. This 
indicated a need for further study of the calculation procedure to determine its general 
applicability. 

Inviscid flow and boundary layer calculations were combined iteratively to attempt to 
account for real flow effects. An existing method was extended to apply to cases with 
boundary layer separation. A second method was investigated which entailed adding a 
calculated boundary layer displacement thickness to the original blade profile to obtain a 
new effective blade profile before recalculating the inviscid flow. The former method 
was judged more accurate and less cumbersome to apply. 



INTRODUCTION 


The real flow through a compressor blade row is compressible, unsteady, three- 
dimensional, highly viscous in certain regions, and may contain both subsonic and su- 
personic flow regions. With the advent of large computers, the ability to analytically 
describe real flow phenomena has been greatly advanced. But it is still impossible, for 
analytical and computer storage reasons, to adequately analyze real compressor flow 
with a single program. However, by making certain simplifying assumptions, such as 
flows that are steady and two-dimensional, with viscous effects confined to a boundary 
layer, the analytical procedure can be divided into separate segments. The results of 
individual calculation procedures can then be combined (e. g. , a two-dimensional inviscid 
flow calculation and a blade surface boundary layer calculation), and approximations to 
the real flow can be made. How closely the computed results approach the real flows 
must be ascertained from comparisons between computed and experimental data. If an 
iterative procedure is required, the degree of accuracy obtained at each stage of the 
computation will determine how many successive applications of the calculation proce- 
dure are necessary, thereby providing a measure of the compromise between accuracy 
and cost. The term cost is used in a qualitative sense to indicate the number of calcula- 
tion cycles and amount of user effort required to effect a solution. The results of com- 
parisons between analysis and experiment will also provide a means for altering and im- 
proving the analytical procedures. 

This report presents a comparison of computed and measured flow parameters for 
a highly loaded compressor stator blade section. The objective is to develop and eval- 
uate several methods for prescribing input boundary conditions to existing analytical 
methods to enable accurate predictions of blade performance to be obtained. The mea- 
ured data were obtained from a conventionally designed solid stator blade (ref. 1). The 
midspan section of the blade, where flow would be expected to approximate two- 
dimensionality, was selected for consideration. All flow velocities were subsonic. The 
data included inlet and outlet station flow properties, blade surface pressures, and total 
pressure loss coefficients obtained from detailed wake surveys. 

The computed results were obtained by using a compressible inviscid flow calcula- 
tion on the blade -to -blade flow plane (ref. 2), a boundary layer method which calculates 
laminar and turbulent compressible boundary layers by integral techniques (ref. 3), and 
a loss calculation which computes loss coefficient using Speidel's method (ref. 4). Of 
particular interest are the calculated loss coefficients and outlet flow angles. Also dis- 
cussed are some problem areas related to the application of these calculation proce- 
dures. 
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ANALYTICAL TECHNIQUES 


This section describes the analytical procedures used in the study and some prob- 
lem areas encountered in applying them. An inviscid flow calculation on the blade-to- 
blade flow plane provides velocities and flow angles. Blade surface velocities from the 
ideal flow calculation are used to compute boundary layer growth. Loss parameters are 
calculated using information from both the inviscid flow and boundary layer calculations. 


Inviscid Flow Calculations 

The equation in the inviscid flow program (ref. 2) which describes the flow in the 
blade -to -blade flow plane is an elliptical partial differential equation for the stream func- 
tion. The inviscid flow program (ref. 2) solves this equation for the subsonic, compres- 
sible flow regime. It is necessary to specify as input the fluid properties, inlet total 
temperature and density, weight flow, blade geometry, inlet and outlet flow angles, and 
the finite difference mesh. As output the program provides blade surface velocities, 
velocities and flow angles at all internal mesh points, and streamline coordinates and 
directions. 

Because the nature of the equations dictates that the solution be of the boundary 
value type, certain variables must be specified on the downstream boundary. In this 
program the outlet flow angle must be specified. This is one of the unknowns of principal 
interest in the compressor design process. Consequently, an additional aspect of the 
inviscid flow problem is to determine a criterion for setting outlet flow angle, probably 
one related to calculated flow conditions in the trailing edge region. 

Historically, inviscid flow programs were first applied to wing sections which had 
sharp or cusped trailing edges. The downstream boundary condition was specified such 
that the stagnation point was located on the sharp trailing edge in order to avoid the 
mathematical condition of infinite velocities around the trailing edge. Differences be- 
tween calculated and measured flow behavior were ascribed to viscous effects. 

Compressor blade profiles, which generally have rounded trailing edges, cannot be 
analyzed in the same manner because there is no way to know the location of the stagna- 
tion point. One method of specifying the downstream boundary condition that has been 
used is analogous to the original wing methods. With this method an attempt is made to 
minimize the velocity differences between suction and pressure surface at the trailing 
edge (analogous to avoiding the infinite velocity condition). This method was used ini- 
tially in the present study. To apply it, the inviscid flow calculation is performed sev- 
eral times, each for a different value of outlet angle. The solution is chosen for which 
the pressures (or velocities) on the suction and pressure surfaces are equal at the blade 
trailing edge. Normally, several iterations are required to find the exit angle that satis - 
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fies the closure requirement. The flow condition corresponding to closure is treated as 
a base or reference condition. 


Boundary Layer Calculations 

The blade surface velocities from the inviscid flow program are used to calculate 
the boundary layer growth using the program of reference 3. In addition to the surface 
velocities, required input includes upstream flow conditions, fluid properties, blade sur- 
face geometry, and, in some cases, initial boundary layer thicknesses. Among the out- 
put provided by the program are the conventional boundary layer thicknesses, form fac- 
tors, wall friction coefficient, and momentum thickness Reynolds number. 

The program uses integral methods to solve the two-dimensional compressible lam- 
inar and turbulent boundary layer equations in an arbitrary pressure gradient. Cohen 
and Reshotko's method (ref. 5) is used for the laminar boundary layer, transition is pre- 
dicted by the Schlichting-Ulrich-Granville method (ref. 6), and Sasman and Cresci's 
method (ref. 7) is used for the turbulent boundary layer. 

A boundary layer which is initially laminar may proceed through normal transition 
to a turbulent boundary layer, or it may undergo some form of laminar separation before 
becoming turbulent. To provide flexibility for analyzing this behavior, several program 
options (discussed subsequently) are available to the user, some of which were used in 
the present study. The calculations may proceed from a laminar boundary layer through 
transition to turbulent calculations. However, if laminar separation is predicted before 
transition, the turbulent calculations may be started by specifying a factor (C^) by which 
the last calculated value of momentum thickness is multiplied. This new momentum 
thickness and a value for form factor based on the last calculated momentum thickness 
Reynolds number are used as initial values for the turbulent calculations. Another op- 
tion permits the user to force transition from laminar to turbulent boundary layer at any 
station specified by the user, provided that transition or laminar separation have not al- 
ready occurred before that station. Another option permits the user to begin turbulent 
calculations at any station by specifying the initial values for displacement and momen- 
tum thickness. 

Boundary layer development over compressor blades is a subject that is not yet 
clearly understood. Although the boundary layer on each surface begins as laminar from 
the stagnation point, the extent of laminar flow before either laminar separation or nor- 
mal transition occurs is not readily determined. Some factors which influence laminar 
boundary layer growth and separation on stationary blades are surface pressure gradient, 
turbulence level, Reynolds number, surface roughness, and, to some extent, the degree 
of three-dimensional flow. It is not often that all of these factors are measured in ex- 
perimental compressor tests. Consequently, there is a lack of clear, authoritative in- 
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formation. Furthermore, because most of the attention given to compressor blade 
boundary layer flow is given to the suction surface boundary layer, very little is pres- 
ently known about the development of pressure surface boundary layers. 

The methods used to calculate compressor boundary layer flow have commonly as- 
sumed the boundary layer becomes turbulent at the beginning of an adverse pressure 
gradient, or, because of the turbulent compressor flow environment, it becomes turbu- 
lent very near the blade leading edge (ref. 8). In either case, the initial turbulent 
boundary layer conditions must be specified, and this information is not known with cer- 
tainty. Although the assumption of early transition to turbulent boundary layer may be 
correct for many, or even most compressor applications, experimental evidence has 
recently been presented by Walker (ref. 9) which shows that laminar boundary layers 
can persist for long distances on actual compressor stator blading (as distinguished from 
cascade tests). 

Consequently, in the present study, boundary layers were calculated as being ini- 
tially laminar from the leading edge region on both blade surfaces. A discussion of the 
results, and of the subsequent steps taken, are presented in the RESULTS AND DISCUS- 
SION section. 

Once a turbulent boundary layer has been established, a criterion must be estab- 
lished for indicating turbulent separation. Sandborn and Liu (ref, 10) have referred to 
turbulent boundary layer separation as "one of the most important unsolved problems in 
fluid mechanics, " and the Stanford Conference Evaluation Committee (ref. 11) concluded 
that no boundary layer calculation method then in use (1968) was mathematically or phys- 
ically legitimate in the separation region, and that data against which to compare analyt- 
ical methods were equally suspect. It is therefore clear that any turbulent separation 
criterion will be approximate and the results must be interpreted accordingly. 

A separation criterion common to compressor blade analyses which use integral 
boundary layer methods is the incompressible form factor H^. Values of 1. 8 to 2. 6 
have been proposed and used in the past (e. g. , von Doenhoff and Tetervin, ref. 12). Ex- 
perience by the author with the boundary layer method of reference 3 on several blade 
shapes has shown a trend for rapid calculated boundary layer growth once has 
reached 2. 0. This indicates that choices of greater than 2. 0 will not appreciably 
alter the location of separation. In addition, it is likely that the equations lose validity 
in regions of high corresponding to separated flow. Therefore, a value of 1^ = 2. 0 
was selected to signal turbulent boundary layer separation. This criterion will provide, 
at best, a consistent engineering estimate of the location of turbulent separation and an 
approximation of the corresponding boundary layer thicknesses at that location. Although 
the experimental data are limited, an attempt will be made in the RESULTS AND DIS- 
CUSSION section to assess the accuracy of this criterion. 
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Loss Calculations 


Various investigators have developed models for calculating total pressure loss 
coefficient based on boundary layer parameters and flow properties at the trailing edge 
of a blade. Lieblein’s incompressible flow analysis (ref. 13) and Stewart’s compressible 
analysis (ref. 14) were accomplished for the case of a boundary layer which remained 
stable to the trailing edge. If boundary layer separation was indicated, an estimate of 
the loss could be obtained only by extrapolating the respective boundary layer thicknesses 
to the trailing edge station. Achieving good agreement with experimentally measured 
values of loss coefficient using this procedure would seem quite unlikely. 

Because the calculation of turbulent boundary layer separation is not exact, the 
prospects for calculating loss in a separated region are not encouraging. For the pres- 
ent, such ventures must be regarded as approximate. Speidel (ref. 4) performed a one- 
dimensional analysis and obtained an expression for the loss due to separation. This 
was incorporated into an expression for overall loss coefficient which closely resembles 
Lieblein's expression. Speidel and Schlichting report good agreement of the expression 
with experimental data (refs. 4 and 8). The method was used in the present study. Be- 
cause reference 4 is in German, the method is developed in appendix B. 


Iterative Use of Calculation Procedures 

The goal of approximating real flow by separately calculating the ideal flow and the 
boundary layer development presumably can be achieved through a series of successive 
applications of each calculation procedure. For example, the inviscid flow calculation 
will provide the surface velocities for a boundary layer calculation. The boundary layer 
calculation will provide boundary layer parameters and input for a loss coefficient cal- 
culation. Cross-interpretation of the separate calculations may lead to refinements and 
recalculation of individual programs with an aim to converge on a solution which closely 
approximates real flow. 

It is apparent that each individual calculation has an effect on the succeeding one, 
and that several cycles through the series of calculation procedures may be involved. 
This can become expensive, both in computer time and user effort, particularly for 
broad analytical studies in which many blade sections and blade rows are involved. Some 
compromise must, therefore, be affected between accuracy and cost. Consequently, in 
evaluating the effectiveness of the methods employed, consideration will be given not 
only to accuracy but also to the number of calculation cycles and amount of user effort 
required. 
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EXPERIMENTAL BLADE 


The measured data used in this report were obtained from the midsection of a con- 
ventional solid stator used in the baseline stage of reference 1. The design of this sta- 
tor is completely described in reference 15 and the profile geometry is shown in figure 1. 
The stator was designed to be representative of the middle stages of a compressor. A 
cross section of the flow path is shown in figure 2. The lack of flow path curvature is 
evident. A double -circular -arc airfoil shape was used. Design equations from refer- 
ence 16 were used to set camber, incidence, and deviation. A summary of geometrical 


TABLE I. - SUMMARY OF DESIGN PARAMETERS 


50-PERCENT SPAN FROM TIP 


Total chamber, cp^, deg 

54. 35 

Total chord, C T , cm 

5. 96 

Maximum thickness to chord ratio, T max /C T 

0. 09 

Setting angle, y, deg 

15.06 

Solidity, a 

1. 35 

Inlet Mach number, M^ 

0. 57 

Diffusion factor, D 

0. 47 

Incidence angle (to mean line), i, deg 

-2.23 

Deviation angle, 6, deg 

12. 12 

Loss coefficient, w 

0.071 

Inlet velocity, Vj, m/sec 

194. 5 

Inlet axial velocity, V z j, m/sec ' 

149.0 

Inlet relative angle, /3'^, deg 

40.0 

Outlet velocity (axial), V2, Vg z , m/sec 

150. 4 


and aerodynamic design data is given in table I. The measured loss against incidence 
angle curve of the midspan section is shown in figure 3 as an example of the perform- 
ance. 


RESULTS AND DISCUSSION 

In this section the computed and experimentally measured results are compared. 
Measured data represent real flow, and efforts are directed toward determining what 
modifications and/or specific inputs to the analytical methods are necessary to produce 
agreement with measured data. The flow parameters to be compared are the blade sur- 
face pressure coefficients, loss coefficient, and outlet flow angle. 
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The results are presented and discussed in the sequential order in which the com- 
puted solutions must be obtained: inviscid flow, boundary layer, and loss calculations. 
In a final section the foregoing results are evaluated and some iterative procedures are 
synthesized. 


Inviscid Flow Calculations 

Three conditions were chosen for detailed analytical investigation: the design point 
(i = -2. 4°) and conditions having both higher and lower than design incidence angle 
(i = 0. 3° and i = -6. 6°). These conditions are indicated as solid circles in figure 3, 
which shows the operating points at which experimental data were taken. 

The blade surface pressure distributions are shown in figure 4 for both the experi- 
mental and calculated cases. As shown in figure 4, it is assumed for the calculated case 
that pressures on the pressure and suction surfaces become equal at the blade trailing 
edge. The calculations compare favorably with experiment over the front portion of the 
blade, but discrepancies become greater toward the trailing edge. In addition, a com- 
parison of the computed and measured exit flow angles (fig. 4) reveals considerable dif- 
ferences. It was found that to close the analytical surface pressure distributions at the 
trailing edge the analysis required 8° to 10° more fluid turning than was measured. 

The differences between measured and calculated results are believed to be due to 
three principal effects. One effect is produced by the uncertainty of the downstream 
boundary condition (placement of stagnation condition on a rounded trailing edge). The 
second effect is that produced by boundary layer blockage and separation. The suction 
surface separation is indicated experimentally for i = -2. 4° and i = 0. 3° by the flat- 
tening of the pressure distribution near the trailing edge region. The third effect is the 
degree of three-dimensional flow actually present in the measured data. At the 
50-percent span location there should be a minimum amount of radial flow. There will 
be some radial flow present, however, and it will affect the comparison with calculated 
flows. It is not possible to quantify any of the three aforementioned effects. 

The inviscid flow calculations were repeated using as input the measured values of 
fluid exit angle and observing the associated trailing edge flow condition. The calculated 
surface pressure coefficients are shown in figure 5 and two general results should be 
noted. First, using the experimental exit angle in the calculations results in a closer 
comparison between experimental and calculated surface pressures, particularly from 
the 50-percent position rearward. Second, using the experimental exit angle resulted 
in the crossing of the computed surface pressure distributions upstream of the trailing 
edge. In fact, the distributions for all operating conditions crossed between 91. 5 and 
92. 5 percent of axial chord. This suggests the possibility that a criterion for prescrib- 
ing exit angle may be devised which is a function of the crossing point (among other pa- 
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rameters) of the surface pressure distributions. Methods for prescribing exit angle 
will be discussed in more detail later in this section. 


Boundary Layer and Loss Calculations 

No experimental boundary layer measurements were taken, thereby preventing di- 
rect confirmation of analytical boundary layer computations. However, exp erimental 
wake surveys were conducted. These permit comparisons of experimental and calculated 
loss coefficients, but they do not permit the source of any discrepancies to be identified. 

The boundary layers were calculated as laminar on each blade surface beginning at 
the stagnation point. In every case, due to the magnitude of the pressure gradient, lam- 
inar separation was indicated before transition on both the suction and pressure surfaces. 
The calculations were then restarted from the indicated point of separation. To begin 
the turbulent boundary layer calculations the initial conditions used were the last calcu- 
lated laminar values for momentum thickness and Reynolds number based on momentum 
thickness. A sample plot of incompressible form factor and displacement and momen- 
tum thicknesses for the design point condition is shown in figure 6. 

The calculated loss coefficient corresponding to the boundary layer distribution 
shown in figure 6 was w = 0. 017. Loss calculations made at off-design conditions, for 
similar treatment of the boundary layer, showed little variation of to with incidence 
angle. The calculated loss coefficients are considerably lower than the experimental 
values of to = 0. 054 to 0. 093 shown in figure 3. 

An additional boundary layer calculation was made in which the boundary layer was 
initially laminar, and after laminar separation occurred the initial turbulent momentum 
thickness was a factor of five greater than the last calculated laminar momentum thick- 
ness before laminar separation. An increase in momentum thickness of this magnitude 
appears to be out of the range of reported boundary layer experience. Despite this, 
however, its effect on calculated loss coefficient was only to increase it to to = 0. 037. 

Because the calculations for initially laminar boundary layers result in very large 
increases in boundary layer thickness through the separation-transition region, an al- 
ternative approach appears necessary. One such approach, as described earlier, is to 
assume the boundary layer to be turbulent from very near the leading edge on both blade 
surfaces as Schlichting (ref. 8) and others have done. This procedure was followed 
next. Some initial considerations must first be discussed, however. 

When the boundary layer is assumed initially turbulent, what is implied is that the 
laminar boundary layer separates and reattaches almost immediately. The important 
unknowns in this problem then include initial turbulent displacement and momentum 
thicknesses. They must be determined by inference - that is, what values of initial dis- 
placement and momentum thickness produce agreement between measured and computed 
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loss coefficient? (Loss coefficient is a direct function of momentum thickness at the 
trailing edge, and, therefore, related to the initial value chosen at the leading edge, 
see eq. (B2). ) This procedure implicitly assumes that the boundary layer and loss 
models are both valid, an assumption that cannot be proved at this time. 

A further consideration which will affect the loss calculations is the relation be- 
tween initial thicknesses chosen for the suction surface and those chosen for the pres- 
sure surface. Because the initial conditions occur in different pressure gradients on 
each surface, and because the gradients change with incidence angle, it is conceivable 
that different initial turbulent boundary layer thicknesses should be chosen for each sur- 
face. However, the testing of the large number of possible combinations, none of which 
are necessarily unique, would complicate and lengthen the computation process. There- 
fore, for the sake of simplicity the same initial displacement and momentum thickness 
used on the suction surface was used also on the pressure surface for any given calcula- 
tion. 

In choosing initial values, 6* and 9 were always chosen such that the initial ratio 
of 6 */9 was 1. 4. There are several precedents for this choice. Schlicting (ref. 6) 
reports that transition on a flat plate results in an initial turbulent form factor of 1. 4. 
Gostelow (ref. 17) used an initial form factor of 1. 4 in his analytical calculations for 
compressor blades in cascade. And McNally (ref. 3) found that on airfoil tests (ref. 18), 
natural transition was calculated, and the initial turbulent form factor was 1. 4. 

One final question relates to the effect of initial thickness on turbulent boundary 
layer separation prediction. Using the Speidel loss model, loss calculations depend on 
separation location (local surface velocity and blade profile thickness at the separation 
point are used) and on boundary layer parameters at separation. If different initial 
thicknesses were found to produce different separation locations, then it might be pos- 
sible to find two or more initial thickness conditions which correspond to the same loss. 
To investigate this effect, a series of suction surface boundary layer calculations were 
made for the design incidence operating condition. A plot of incompressible form fac- 
tor is given in figure 7 for four initial thicknesses ranging from 0. 00609 to 0. 0609 cen- 
timeter. Differences are seen over the early portion of the blade where form factor is 
not critical. But all curves show the same behavior over the latter portion of the blade 
where form factor rises sharply to separation values. The location of indicated separa- 
tion shows little variation for the four initial thicknesses considered. 

Analytical calculations were made for the three operating conditions previously 
noted (i = -6. 6°, -2. 4°, and 0. 3°). One additional operating point was considered for 
this set of calculations, the condition where i = 2. 88° (see fig. 3). Boundary layer and 
loss calculations were performed for a range of initial thickness values until loss coef- 
ficients were brought into agreement with experimental values. The initial displacement 
and momentum thicknesses which produced this agreement at each incidence angle are 
given in table II. Plots of the calculated turbulent boundary layer distributions at each 
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TABLE II. - INITIAL BOUNDARY LAYER THICKNESSES AND 


CORRESPONDING ANALYTICAL LOSS COEFFICIENTS 


Incidence 

angle, 

i, 

deg 

Initial 

displacement 

thickness, 

S’ 

cm 

Initial 

momentum 

thickness, 

5 i> 

cm 

Analytical 

loss 

coefficient, 

u> 

Experimental 
loss coeffi- 
cient, 
w 


0.0914 

0. 0652 


BM 


.0685 

.0490 




. 0610 

.0435 

. 064 

. 062 


.0426 

.0345 

. 056 

.055 


incidence angle are presented in figure 8. Separation is indicated when the incompres- 
sible form factor curve intersects the line = 2. 0. 

Two results should be noted. First, in order to predict the experimental loss coef- 
ficient at the design point, i = -2. 4, it was necessary to begin the turbulent boundary 
layer calculations with an initial 6* of 0. 061 centimeter and an initial 6 of 0. 0435 cen- 
timeter. These are rather large values, corresponding to 1 and 0. 7 percent of blade 
chord, respectively. The second point to be noted is that the initial boundary layer 
thicknesses increased as incidence angle increased. 


Prediction of Turbulent Boundary Layer Separation 

Since no experimental boundary layer measurements were taken, it was therefore 
not possible to establish the location of boundary layer separation accurately. The only 
measurements available are surface pressure distributions. The criterion that was 
used to locate separation was to note when this distribution changed slope and became 
flat.' It was assumed that when separation occurred velocity deceleration ceased and 
static pressure remained constant. This will provide a gross estimate of the location 
only, since measurements to indicate streaking, velocity profile reversals, and shear 
stress distributions are not available. 

To illustrate the process, refer to figure 4(a) in which a dashed line at constant Cp 
is drawn through the last experimental point on the suction surface. At the X/C^. loca- 
tion where this line intersects the faired curve through the nonseparated points (dashed 
line also), separation was inferred to occur. 

This procedure is clear for the conditions shown in figures 4(a) and (b). However, 
in figure 4(c) (i = -6.6) suction surface separation is not clearly indicated. It was in- 
ferred to occur at the last experimental point (X/C^ = 0. 86) by using the following process. 
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From the experimental data presented in figures 4(a) and (b), the difference in suc- 
tion surface Cp from the inferred separation point to the minimum value is as follows: 

(a) For i = 0. 3, AC p = 0. 17 - (-0. 5) - 0. 67 

(b) For i = -2. 4, AC p = 0. 13 - (-0. 53) = 0. 66 

In each case, this corresponds to a maximum surface velocity to velocity at separation 

of V /V = 1. 33. Although the last suction surface data point in figure 4(c) does 
s op 

not show evidence of flattening, the AC p at that point is 0. 69 and V max /V local = 1. 33. 
Therefore, it was felt that there was adequate reason to use it as an indicated separation 
point. 

Taking the axial chord locations thus determined, and going to the analytical bound- 
ary layer curves of H. against percent axial chord (fig. 8), one can find the corre- 
sponding IT values at separation. These are shown in table III. They are relatively 
close to the = 2. 0 values used in this study. Although this is a rather unrefined 
process, it does provide some justification for the use of H i = 2 . Q as a separation cri- 
terion. Additional data and more refined methods of evaluating experimental separation 
location will aid in determining the general applicability of this criterion. 


TABLE III. - FORM FACTORS AT SEPARATION 


Incidence 

angle. 

i, 

deg 

Pressure 

coefficient 

rise. 

Axial loca- 
tion of 
separation. 

x/c t 

H. at 
separation 
(from fig. 7) 

0. 3 

0.67 

0. 765 

1. 88 

-2. 4 

.66 

.815 

1. 88 

-6.6 

. 69 

. 86 

1. 94 


Iterative Applications of Inviseid and Viscous Calculations 


In the foregoing sections the results obtained from once-through applications of the 
calculation procedures have been presented. In order to approximate real flow per- 
formance some refinements to the procedures are clearly necessary, and indeed, were 
anticipated. In this section some approaches involving iterative applications of the cal- 
culation procedures will be developed and discussed. 

General observations . - Inviseid flow calculations presented thus far allow several 
conclusions. First, choosing an exit flow angle such that surface pressure distributions 
closed at the trailing edge predicted, as expected, much more fluid turning that was 
actually achieved. Also, surface pressures began to diverge from measured values as 
the trailing edge region was approached. Second, using the experimentally measured 
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exit angle in the analytical calculations resulted in good agreement between experimental 
and analytical surface pressures. It is instructive to know that if the exit flow angle can 
be approximated, the surface pressures resulting from a single pass through the invis- 
cid flow calculations are in reasonable agreement with measured values. The search 
then becomes centered on methods for approximating exit flow angle. In this regard, it 
is of interest to note that in the calculations using measured exit angle (fig. 5) all of the 
analytical surface pressure distributions showed crossing between 91. 5 and 92. 5 percent 
of axial chord. 

This therefore suggests that one possible approach may be to construct a relation 
between the position of crossing of the surface pressure distributions (and thereby, the 
exit angle that produced it) and the pressure gradient or other aerodynamic parameters 
which are available from the inviscid flow calculations. The correlation with aerody- 
namic parameters available from the inviscid flow calculations would be desirable in 
order to avoid complex empirical relations relating blade shape, thickness distribution, 
solidity, camber, and other geometric parameters. If such a relation could be con- 
structed, it would permit the exit angle to be determined from application only of the 
inviscid flow calculation thereby reducing analytical complexity and cost. However, the 
synthesis of such a relation requires additional data from blades of widely differing de- 
signs. 

Gostelow method extended . - Gostelow (ref. 17) was successful in approximating 
real flow by using only an inviscid flow calculation. However, he investigated relatively 
lightly loaded blades (C4 profile with 30° camber) in low subsonic flow. He predicted 
exit flow angle and surface pressure distributions quite closely using the Martensen in- 
viscid flow method and a "fairing-in” process near the trailing edge of the blade. 

The "fairing-in” process is an empirical method and constitutes an approximation 
for viscous and rounded trailing-edge effects. Several inviscid flow calculations are 
made, and a tangent to the C p distributions is faired in at the 85-percent chord loca- 
tion on each surface (see fig. 9). The "correct" solution (exit angle condition) is the 
one for which the faired-in tangents just meet at the trailing edge (100-percent chord). 
Gostelow followed Preston’s analysis (ref. 19) which asserts that no net vorticity is dis- 
charged into the wake and that wake curvature effects may be neglected. The 85-percent 
chord location was chosen on the basis of work reported by Spence and Beasley (ref. 20). 

It should also be noted that, in most cases, for the condition corresponding to the 
faired-in tangents meeting at the trailing edge, the calculated surface pressure curves 
will cross inside the trailing edge. This agrees with the results obtained from calcula- 
tions in this study, as noted in the previous section. 

In the examples shown in reference 17 the suction surface boundary layer was well 
behaved, showing no evidence of separation. It seems unlikely that the method can be 
applied directly to the data evaluated in this report because separation is indicated over 
roughly the last 20 percent of the suction surface. Indeed, an application of the refer- 
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ence 17 method to the plots of figure 5, where measured exit angles were used as input, 
showed the tangents to cross at axial chord locations between 93 and 94 percent rather 
than 100 percent. Thus, a direct application of Gostelow’s method would have overes- 
timated the fluid turning in these cases. 

A method has been devised for modifying this general approach to permit applica- 
tion to blades having separated boundary layers. In brief, the method applies the tan- 
gent criterion at the indicated separation point when it occurs upstream of the 85-percent 
chord position. More specifically, the method requires iterative use of the inviscid 
flow and boundary layer calculations and is now outlined: 

(1) Perform an inviscid flow calculation iteratively to determine the case for which 
suction and pressure surface Cp curves cross at trailing edge (base condition). 

(2) Calculate the boundary layer assuming it to be turbulent from near the leading 
edge. 

(3) Determine the location of predicted boundary layer separation using H. = 2. 0. 

(4) From the Cp against X/Cj_ distribution for the base condition, determine the 
Cp at the separation location. Subtract the minimum C p from it to get the AC p that 
produced separation. 

(5) Recalculate the inviscid flow at some /3jj > 02base ^ e ' g ' ’ @2 “ ^2base ~ 5 °)- 

(6) Add ACp to the minimum C p to locate the new C p corresponding to separa- 
tion and from which the tangent is to be faired in. Fair in the tangent from this point. 

(7) Fair in the tangent on the pressure surface C p curve from the X/C t = 0. 85 
point. 

(8) Repeat steps (5) to (7) until the faired-in curves cross at the trailing edge 
(X/C t = 1. 0). The exit flow angle corresponding to this final calculation is the approxi- 
mation to the real flow angle. 

The tangent to the suction surface C is faired in from the point of estimated sep- 
aration. If AC is used as the parameter in steps (4) and (6) the pressure which cor- 
responds to separation in the first iteration is transferred to subsequent iterations. Re- 
peated applications of the boundary layer calculation are therefore avoided. The X/Cj. 
was not used as the separation-determining parameter because the point at which separa- 
tion occurs is largely a function of pressure rise (i. e. , ACp) rather than distance. 

The previously discussed compressor stator data are used to check the method. 

Steps (1) to (3) have already been performed (figs. 4 and 8). The series of iterations 
outlined in steps (5) and (8) are not carried out, however. Instead the calculations per- 
formed using experimentally measured exit angle (fig. 5) are used. This calculation 
represents the "correct" condition and if fairing in from the appropriate position results 
in crossing at X/C^. = 1. 0 the method outlined previously is confirmed. 

The boundary layer curves of figure 8 can be used to find the location of suction sur- 
face separation (using IT = 2. 0). From the inviscid flow calculations (fig. 4) at those 
X/Ct locations the Cp at separation and the rise in C p from minimum can be deter - 
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TABLE IV. - DETERMINATION OF 


PRESSURE RISE TO SEPARATION 
FOR FIRST ITERATION 


Incidence 

angle, 

i, 

deg 

Suction surface 

(x/c t ) 

v y sep 

( Cp )sep 

AC 

P 

0. 3 

0. 785 

0. 07 

0. 66 

-2.4 

. 83 

.06 

.66 

-6.6 

. 865 

. 05 

.65 


mined (AC p ). These are given in table IV. 

The ACp values thus determined are applied to the suction surface of the inviscid 
flow calculation of figure 5 as shown in figure 10. Tangents are faired-in from the re- 
spective suction surface points and from the point on the pressure surface corresponding 
to X/Cj. = 0. 85. The faired-in curves intersect at just under the X/C^. = 1. 0 position, 
indicating that the exit angles used in the calculation (which were the experimentally 
measured angles) were very close to what would have been determined by direct applica- 
tion of the modified Gostelow method. The locations of the intersections of the faired-in 
tangents as determined from figure 10 are as follows: 


Incidence 

Location of 

angle, 

intersection 

i. 

of tangents, 

deg 

x/c t 

0. 3 

O'. 99 

-2. 4 

. 99 

-6.6 

. 98 


Therefore, the modified Gostelow method produced an accurate result, at least for 
the set of data evaluated. This approach involves several iterations of the inviscid flow 
calculation and a single pass through the boundary layer calculation. One area of uncer- 
tainty is the choice of initial boundary layer thickness. It was shown earlier (fig. 7) that 
for a typical surface pressure distribution the predicted location of separation was not 
sensitive to initial thickness. The user therefore has some latitude in the choice of 
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initial 6*. Once 6* is specified, an initial form factor of 1. 4 is appropriate and will 
then specify momentum thickness. 

Displacement-thickness-added method . - The method just described approximates 
the effect of the boundary layer by altering the inviscid flow surface pressure distribu- 
tion. Another method for approximating the effect of boundary layer is to add calculated 
boundary layer displacement thickness to the blade profile geometry and recalculate the 
inviscid flow over the redefined blade. The process can be repeated until convergence 
within desired limits is achieved. An outline of the procedure is as follows: 

(1) Perform the inviscid flow calculation iteratively to determine the case for which 
suction and pressure surface Cp curves cross at the trailing edge (base condition or 
ideal condition). 

(2) Calculate the boundary layer, assuming it to be turbulent from the leading edge. 

(3) Add the boundary layer displacement thickness to the original blade profile, thus 
defining a new effective blade profile. 

(4) Recalculate the inviscid flow for the new effective profile. 

(5) Repeat steps (2) to (4) until inviscid flow results converge within acceptable 

limits. I \ 

This procedure was carried out for the design point (i = -2. 4). The inviscid flow 
calculation shown in figure 4(b) constitutes the base condition (step (1)). The boundary 
layer calculation corresponding to it, shown in figure 8(b) (initial displacement thickness 
on both surfaces, 0. 061 cm), constitutes step (2). The calculated displacement thick- 
ness was added to the blade profile geometry. Since boundary layer separation was in- 
dicated on the suction surface, the new profile was extrapolated from the separation 
point to the trailing edge to continue a smooth profile shape. The inviscid flow calcula- 
tion requires geometry with a trailing edge radius, so the profile was closed with a cir- 
cular arc. The resulting blade section profile is shown in figure 11. 

The pressure distribution calculated by the inviscid flow program is shown in fig- 
ure 12. The local accelerations near the trailing edge are caused by the large trailing 
edge radius. Because of the local accelerations it was necessary to devise a method for 
setting exit angle which was not dependent on the calculated surface pressures. In the 
approach used, it was assumed that the flow follows the direction of the effective profile 
(original geometry plus boundary layer displacement thickness) to the point at which the 
trailing edge circle joins the profile. The downstream angle is then taken to be an aver- 
age of the two surface angles at the respective points of tangency with the trailing edge 
circle. This procedure neglects three-dimensional effects and assumes complete mix- 
ing to a uniform downstream angle. If separation is predicted (as it was in this case), 
fluid turning is assumed to cease at the predicted separation point and the surface angle 
there is used in the average. The surface angles of the new effective blade are directly 
available from the spline fits in inviscid flow calculations. 
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It is clear from the results that because of the large boundary layer thickness and 
separation the method did not converge in this case. In this, and similar cases, one 
pass through the calculation procedure is probably sufficient. 

In general, this method should converge. The reason it did not in this case is that 
the initial boundary layer thickness was chosen in order to match calculated loss with 
measured loss. 

The exit angle obtained by this procedure was -1. 25° compared to 1. 6° measured 
experimentally. The surface pressure distributions obtained (fig. 12) show very good 
agreement with surface pressure data on the pressure surface. Although the predicted 
pressures on the suction surface are displaced from the data somewhat, agreement is 
acceptable. 

Exit condition adjusted for separation . - One additional iterative method has been 
conceived but has yet to be evaluated. The data from figure 4(a) are reproduced in fig- 
ure 13 and dashed lines have been faired through the points. Boundary layer separation 
on the suction surface is indicated by a flattening of the pressure distribution near the 
trailing edge and is drawn as a constant C p line. It can be observed that when the 
separated region is represented as a constant pressure region, the suction and pressure 
surface distribution appear to close at the trailing edge (100-percent chord location). 

That is, the pressure surface flow adjusts to the region of higher static pressure pro- 
duced by the suction surface separation near the trailing edge. 

A pressure distribution similar to figure 13 could be produced analytically and would 
thereby provide a means of establishing exit angle. The procedure would require the 
following steps. Input surface pressures for a boundary layer calculation are provided 
by an inviscid flow calculation. The location (X/C t ) of suction surface boundary layer 
separation is predicted by the boundary layer calculation. The value of on the suc- 
tion surface at the separation location (X/C J is extended at constant C to the trailing 
edge. The pressure surface Cp distribution is extrapolated to the trailing edge. If the 
suction surface and pressure surface Cp values are equal at the trailing edge station, 
the /?2 used in the inviscid flow calculation is the predicted value. If not, the series of 
calculations is repeated until pressure surface and adjusted suction surface pressures 
agree at the 100-percent chord location. 

Comparison of analytical methods . - Several approaches have been suggested for 
utilizing inviscid flow and boundary layer calculation procedures to approximate real 
flow conditions. Particular attention has been paid to the ability to predict surface pres- 
sure distributions and outlet angle. In this section some advantages and disadvantages 
of each approach are enumerated to allow a comparative assessment to be made of the 
utility (cost) and accuracy of the suggested approaches. 

One approach, yet undeveloped, would involve a correlation of the location of cross- 
ing of the suction and pressure surface pressures (and thereby the exit angle that pro- 
duced it) and aerodynamic parameters available from the inviscid flow calculation and 
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perhaps the boundary layer calculation as well. Such a method appears feasible. It 
would be very simple to apply. But much data correlation would be required to develop 
it. 

The extended Gostelow method showed a good ability to predict exit angle. Predic- 
tions of surface pressure distributions were good, but not exact. (In fig. 12 the extended 
Gostelow method and the displacement thickness added methods are compared, and their 
abilities to predict surface pressure distribution appear to be about equal. ) The basic 
method developed by Gostelow is well documented and was shown to work well for low 
cambered blades (refs. 17 and 21). The extension in this report to highly cambered 
blades showed equally good promise. It requires several iterations to apply but consists 
of relatively simple operations. It retains the disadvantages of the basic method, that 
is, the empirical basis, and the necessity to manually fair-in tangents. 

The displacement-thickness -added method predicted surface pressures with about 
equal accuracy as the extended Gostelow method. Predictions of fluid exit angle were 
acceptable. The method at present does possess some serious disadvantages, however. 
The assumptions for initial boundary layer thicknesses in the boundary layer calculations 
are crucial in this method. The initial thicknesses will determine the new profile geom- 
etry: profile thickness, local surface slopes, and trailing edge thickness. And there is 
little reliable information to guide the user in this choice. Furthermore, the necessity 
to fit a spline curve through the new effective profile coordinates is time consuming, and 
it must be done correctly to obtain meaningful results. 

Weighing the advantages and disadvantages of each method, based on the experience 
gained in this study, the extended Gostelow method appears to be the more serviceable 
method. It demonstrates acceptable accuracy and, although several iterations are nec- 
essary to arrive at a solution, the procedures are not cumbersome. 

The two iterative methods discussed showed reasonably close comparison between 
calculated and measured blade surface pressures (velocities). The user is therefore 
able to obtain a direct and reasonably accurate determination of blade surface velocity 
deceleration. It is interesting to note that the deceleration of blade surface velocity has 
previously been utilized to give a measure of the blade loading. Lieblein's D-factor 
(ref. 21), for example, is an attempt to calculate the velocity deceleration on the blade 
suction surface from the known velocity diagrams at the blade inlet and outlet stations. 
This D-factor was developed for near design incidence angle applications. The ability 
to calculate blade surface velocities (pressures) with reasonable accuracy offers a more 
direct method to indicate blade loading levels at all operating conditions and, therefore, 
to predict useful operating range. 
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SUMMARY OF RESULTS 


Analytically computed flow parameters for the midspan, double-circular arc blade 
section of a highly loaded compressor stator have been compared with measured flow 
parameters. Analytical procedures used included calculations for inviscid flow on a 
blade -to-b lade flow plane, for boundary layer on the blade suction and pressure surfaces, 
and for total pressure loss coefficient. All velocities were subsonic. Comparisons 
were made at three incidence angles between experimentally measured and computed 
blade surface pressure distributions, outlet flow angle, and loss coefficient. 

Two-dimensional calculation methods proved adequate for predicting performance 
at the midspan location where three-dimensional effects are believed to be minimized. 
Radial flows cannot be considered completely absent, however, and some portion of the 
discrepancy between calculated and measured flows must be attributed to them. 


Inviscid Flow Calculations 

In order to deduce the relations between outlet flow angle and calculated flow condi- 
tions at the blade trailing edge region, the inviscid flow calculations were carried out by 
two methods. In the first, the pressures on the blade suction and pressure surfaces 
were made equal at the blade trailing edge station and the associated fluid outlet angle 
accepted as a result. In the second; the measured outlet flow angles were used as input 
and the associated calculated flow conditions in the blade trailing edge region accepted 
as a result. Comparison of calculated and measured data gave the following results: 

1. When the exit flow angle was adjusted to make calculated suction and pressure 
surfaces pressures equal at the blade section trailing edge, the following observations 
were made: 

(a) Favorable agreement between analytical and experimental pressures was 
achieved over the front portions of the blade, but divergence increased toward the 
trailing edge where the boundary layer and separation effects are greater. 

(b) The analytical outlet angles differed significantly (as much as 10° greater 
turning) from measured values. 

2. When the experimentally measured outlet angle was used in the inviscid flow cal- 
culations, much better agreement between experimental and analytical surface pressures 
was achieved, although some differences still existed in the trailing edge region. For 
all three incidence angle conditions the calculated suction and pressure surface pres- 
sures became equal at the 91. 5 to 92. 5 percent chord location (i. e. , surface pressure 
distributions crossed inside the trailing edge). 
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Boundary Layer and Loss Calculations 


Boundary layer and loss calculations yielded the following results: 

1. Laminar boundary layer calculations predicted laminar separation to occur 
before transition on both blade surfaces at all incidence angles. 

2. Transition to a turbulent boundary layer was forced analytically at the indicated 
laminar separation point, but agreement between analytical and measured losses was 
not achieved, even when the boundary layer thickness was allowed to increase by a fac- 
tor of five. 

3. All boundary layers were subsequently calculated as initially turbulent from near 
the blade leading edge. Calculated loss values were made from a flow model developed 
by Speidel which accounts for boundary layer losses to the separation point and includes 
an additional term to account for loss due to separation. To make calculated losses 
match measured losses the following had to be done: 

(a) At the design point, initial turbulent boundary layer displacement thickness 
had to be 0. 061 centimeter and initial momentum thickness had to be 0. 0435 centi- 
meter on each blade surface. The large required values of initial thicknesses indi- 
cate the need for further application of the calculation procedure to determine its 
general applicability. 

(b) As incidence angle increased, the initial thickness values had to increase. 


Iterative Calculations 

Two methods were considered which involved iterative use of inviscid flow and 
boundary layer calculations. Each produced acceptable results: 

1. The first method introduced a modification (accounting for separated boundary 
layers) to a procedure reported by Gostelow. Boundary layer and rounded trailing edge 
effects are approximated by neglecting calculated surface pressures near the trailing 
edge and fairing-in tangents to the pressure distributions at a location upstream of the 
trailing edge such that the tangents cross at the trailing edge station. By iterative ap- 
plication of this method fluid exit angle was determined quite closely and surface pres- 
sures showed good agreement with experimental data. The boundary layer calculations 
were used to locate the point of turbulent separation from which a tangent was faired in 
to the pressure distribution. 

2. The second method investigated entailed adding calculated boundary layer dis- 
placement thickness to the original blade profile geometry to define a new "effective" 
blade shape and then recalculating inviscid flow around this new blade profile. Resulting 
surface pressures showed good agreement with measured data. The prediction of the 
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fluid exit angle was less accurate than the extended Gostelow method, and in general, 
the method was more cumbersome to apply. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, June 29, 1973, 

501-24. 
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APPENDIX A 


SYMBOLS 


blade chord length, cm 


static pressure coefficient, 


diffusion factor, 1 



V 


P l “Pi 



91 Y 9 2 
2aV' 1 


force, N 

incompressible form factor (ref. 3) 

incidence angle, i. e. , angle between inlet flow direction and line tangent to blade 
section meanline at leading edge, deg 

Mach number 

o 

static pressure, N/m 
blade spacing, cm 
blade thickness, cm 
velocity, m/sec 

distance from blade meanline to blade surface (fig. 14), cm 
axial distance, cm 

fluid angle, with respect to axial direction, deg 
fluid turning angle, /3k - i3p deg 

deviation angle, i. e. , angle between exit flow direction and line tangent to blade 
section meanline at trailing edge, deg 

boundary layer displacement thickness, cm 

blade setting angle, deg 

boundary layer momentum thickness, cm. 

3 

fluid density, kg/m 

blade solidity, chord to spacing ratio 

blade camber angle, deg 

blade loss coefficient, ratio of total pressure loss to inlet dynamic pressure 



Subscripts: 
l local 

max maximum 

p pressure 

s suction 

T total 

z axial direction 

9 tangential direction 

1 inlet flow plane 

2 outlet flow plane 
Superscript: 


relative 



APPENDIX B 


DEVELOPMENT OF SPEIDEL’S LOSS MODEL 


Speidel (ref. 4) followed a similar line of development as Lieblein (ref. 9) to arrive 
at a loss coefficient expression for unseparated flow expressed in terms of boundary 
layer parameters. The analysis was for incompressible flow. The resulting expression 
for loss coefficient (in the nomenclature of the present report) was 


a 


(Bl) 


where 0 C is the momentum thickness on the suction surface at the trailing edge and 9 
is the corresponding thickness on the pressure surface. 

When separation occurs, presumably on the suction surface, Speidel expressed the 
additional loss in terms of a hypothetical momentum thickness 6^. The rationale used 
to achieve an expression for 9 ^ follows. 

Figure 14 represents the blade thickness against distance near the trailing edge re- 
gion. The thickness distribution shown is symmetrical around the mean line. Point A 
on the suction surface is the point of separation and the distances from the mean line to 

the suction and pressure surfaces at that point are y 0 and y . Plotted below the 

S A PA 

thickness distribution plot is a representation of surface pressure against distance. The 
nonseparated flow case is shown in solid lines. Both surface pressures meet at the 
trailing edge at pressure Pjj. At point A the flow separates and the static pressure re- 
mains constant to the trailing edge. It is assumed that the pressure on the pressure 
surface accommodates itself to the new trailing edge pressure p A in a linear fashion. 

It is also assumed that main flow total pressure remains constant between A and H. 

Shown in figure 15 is a plot of surface static pressure against blade thickness. The 
nonseparated flow case is once again shown with solid lines and the separated case in- 
dicated by dotted lines. In the latter case, it can be seen that as a condition of zero 
thickness (trailing edge) is approached, both pressures become equal to p A . Shown as 
a cross-hatched area is the difference between the separated and nonseparated cases. 
This cross-hatched area corresponds to a defect of force on the blade due to the bound- 
ary layer separation, and it can be expressed as a loss of momentum. 

The defect in force can be expressed as the area of the cross-hatched triangles in 
figure 15. On the suction surface of the mean line, 
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and on the pressure surface side, 


F s ■ <Ph ■ Pa> \ ?S A 


F P “ ( p H ' P A* 2 y p A 


The total force defect due to pressure is 

F T = Fg + Fp = (p H - P A )y A 


where 

y A ' y S A + y P A 

and 




r %-» 

p A 


The incompressible Bernoulli equation expressed between A and H, with total pressure 
assumed constant, is 


P TOT " P A + 2 W A “ P H 



P H 



so that now 


F T=f W i- W H y A 

The defect in momentum created by separation can be expressed in terms of a mo- 
mentum thickness 0 A : 


2 

Defect in momentum = 
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where W 2 is the average velocity in the outlet plane. The defect in momentum is 
equivalent to the total force defect, that is, 


F T = P*A W 2 


Substituting for F T gives 


? K - w l 


)y A - 


A w 2 


A 2 


,W a\ 2 /V* 


L\ w 2 


W r 


Speidel reports that experimental measurements in the region of the trailing edge show 


so that 



The expression for loss coefficient when separation has occurred is then 


- _ 2 2 ( 0 S + 0 p + 9 d / cos 0*1 
s cos $2 Icos @2 

where 9 is now the momentum thickness at the point of separation and is the mo- 
mentum thickness at the trailing edge. 
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Meridional distance, M, cm 

Figure 1. - Stator blade section geometry at midspan 
section. 



Figure 2. - Cross section of flow path showing lack of curvature. 


• 12 1 


13 

c? .081— 


.04K 


• Flow conditions at which detailed 
analytical calculations performed 

O o 


O 








0L 

-15 


-10 -5 0 

Incidence angle, i, deg 


■ Figure 3. - Stator experimental loss coefficient at mid- 
span section. 
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Axial distance, z/c^ z 

(c) Incidence angle, i - -6. 6°. 

Figure 4, - Comparison of experimental and analytical surface pressure coefficients. 
Closure at trailing edge. 
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(c) Incidence angle, i ■ -6.6°. 


Figure 5. - Comparison of experimental and analytical surface pressure coefficients. 
Experimental exit angle used in analytical calculations. 




Figure 6. - Boundary layer development with initially laminar boundary. Incidence angle, i ■ -2.405. 







Suction surface incompressible form factor, Hj Suction surface mo- 

mentum thickness, 



Axial distance, z/cj z 

Figure 7. - Effect of initial turbulent boundary layer thickness 
on incompressible form factor and momentum thickness. 
Analytical pressure distribution; design incidence, i * 

-2. 405°. 
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Figure 9. - Illustration of "fairing-in" prxess (from Gostelow, 
ref. 17). 



(a) Incidence angle, i =0.3°; exit fluid angle, =1.6°. 



Figure 10. - Extended Gostelow method using faired-in tangents in trailing-edge region. 



Original blade section 



Figure 11. - Stator blade section with boundary layer 
displacement thickness added. 



Tangents 


Figure 12. - Comparison of experimental and analytical surface 
pressure coefficients. Corrected and uncorrected for bound- 
ary layer design point, i - -2.4°. 
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Figure 13. - Experimental surface pressure coefficients 
showing effect of boundary layer separation on trailing 
edge closure condition. Incidence angle, i = 0. 3°. 



Figure 14. - Schematic of trailing edge blade 
profile and pressure distribution with and 
without boundary layer separation. 



Figure 15. - Schematic of pressure defect due to boundary 
layer separation. 
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